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Abstract 

We study different models of radiating slowly rotating bodies up to 
the first order in the angular velocity. It is shown that up to this order 
the evolution of the eccentricity is highly model-dependent even for very 
compact objects. 
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1 Introduction 



In recent papers by Gupta et al [0, ^ , it has been shown that for a specific model 
of a slowly rotating object (up to the second order in angular velocity), the 
ellipticity attains a maximum for a surface gravitational potential, ^ ~ 0.18 
(in geometric units). 

This effect discovered by Chandrasekhar and Miller ||, and initially asso- 
ciated with frame dragging, was later related to the reversal of the centrifugal 
force in ultracompact objects. @, |^]. 

In this work we study different models of radiating, (up to first order) slowly 
rotating objects. The motivation to undertake such a task is threefold: 

• firstly, we want to check if the effect reported by Gupta and collaborators 
(ij |^], appears at first order. In this case the reversal of the centrifugal 
force (a second order effect), although may contribute to, has to be ruled 
out as the main cause of the maximum in the eccentricity. 

• secondly, we would like to clarify if there exists a general tendency (model 
independent) in the evolution of eccentricity as the surface gravitational 
potential increases (see also Q for a discussion on this point). 

• finally, it is worth mentioning that the models presented in |Q ^| are only 
valid within the quasi-stationary approximation. In other words, their 
results are obtained by identifying the "history" of the evolution of the 
compact object with a sequence of stationary states. Instead we shall 
consider here systems which depart from hydrostatic equilibrium 

As we shall see below, even for those cases of highly compact matter con- 
figurations, it appears that there is no model independent tendency for the 
eccentricity to attain a maximum. In one set of our simulations models start 
and end in a stationary state. It is obvious that in these cases the eccentricity 
reaches a constant value, but its value depends on the initial data considered. 
Several other simulations were preformed using two different equations of state. 
The models that emerge never reach stationary state. In these cases the eccen- 
tricity grows up continuously, even for surface gravitational potential as large 
as ^ w 0.31 . 

This paper is organized as follows. In the next section we introduce the 
conventions and give a brief description of the method employed to simulate 
the evolution of slowly rotating objects. Section 3 presents the three families 
of models considered. Finally, in the last section our results are discussed and 
contrasted with those obtained by Gupta and collaborators [li} ||. 
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2 The Formalism 



Let us consider a nonstatic and axially symmetric distribution of matter formed 
by a fluid and radiation. The exterior metric, in radiation coordinates ||, takes 
the Kerr-Vaidya form M : 



ds 2 



1 



2mr 



r 2 + a 2 cos 2 9 



+4a sin ( 



sin 



du + 2dudr - 2a sin 8drd4> 



dud(j) — (r 2 + a 2 cos 2 



2 + a 2 cos 2 9 

9 2mra 2 sin 2 9 



(1) 



a 2 cos 2 9 



Here, a is the Kerr parameter, representing angular momentum per unit mass 
in the weak field limit, and m is the total mass. It is worth mentioning at 
this point that the metric above is not a pure radiation solution and may be 
interpreted as such only asymptotically ||. A pure rotating radiation solution 
may be found in reference |)). As we shall show below, although the interpreta- 
tion of the Carmeli-Kaye metric is not completely clear, the model dependence 
of the considered effect is independent of the shape and the intensity of the 
emission pulse, and may be put in evidence even for a tiny radiated energy, 
AM ra j = 10~ 12 Af(0), which for any practical purpose corresponds to the 
Kerr metric. The interior metric is written as [[10[ 



ds 2 = 



e 2p j-^du 2 + 2dudr j - (r 2 + a 2 cos 2 9)d9 2 
+2ae 2(} sin 2 9 j 1 - — j dudcj) - 2e 2f3 a sin 2 9drdcj) 



(2) 



sin 



,V 



T +a 2 + 2a 2 sin 2 9— 
r 



In equations ({j]) and (@), 



null coordinate and 9 



= x° is a time like coordinate, r — x 1 is the 
x 2 and (f> = x 3 are the usual angle coordinates. The 
u-coordinate is the retarded time in flat space-time and, therefore, u-constant 
surfaces are null cones open to the future. 

The Kerr parameter for the interior space-time (|^) is denoted a and, for 
the present work and for sake of simplicity, will be considered constant and 
relevant only (as well as a in eq. (Q)) up to the first order. Notice that in 
these coordinates the r — constant = r s surfaces are not spheres but oblate 
spheroids, whose eccentricity depends upon the Kerr parameter a and is given 

by 

e 2 = 1 - . ( 3 ) 

r 2 + a 2 
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with r s representing the shell where the eccentricity is evaluated. Observe that 
eccentricity as defined by (||) , although it yields the correct Newtonian limit and 
is the natural definition in the context of metrics (Q) and (||), is not invariantly 
defined . However the main argument to use it here, stems from the fact that 
it is that parameter the one whose evolution is studied in references pi ^] . 

The metric elements f3 and V in eq. (|J), are functions of u, r and 6. A 
function rh(u, r, 6) can be defined by 

v = e2 ,f_2n^ 1 ey\ 

\ r 2 +a 2 cos 2 9 J 

which is the generalization, inside the distribution, of the "mass aspect" de- 
fined by Bondi and collaborators and in the static limit coincides with the 
Schwarzschild mass. 

In this work the modeling is performed by means of a general method pre- 
sented a few years ago [ |l2| , which allows the generation of axially symmetric 
slowly rotating (up to the first order), radiating solutions from known static 
"seed" solutions. 

Only a very brief description of the method is given here, we refer the reader 
to Ul for details. 

The method starts by defining two auxiliary functions 

P = TT^- (5) 

and 

P= P - pu ' . (6) 

which are called the effective variables. In equations (||) and (||) lo x represents 
the radial component of the velocity of a fluid element as measured by a local 
minkowskian observer. Also in the above equations the physical (local) pressure 
and density are denoted by p and P. 

Next, it can be easily shown using the field equations that the metric ele- 
ments m and (3 can be expressed as 



m 



/ dr47rr 2 /5 (7) 
Jo 



and 

where r = a(u) is the equation of the boundary surface. 

It is clear that, if the radial dependence of p and P are borrowed from the 
"seed" static solution, the metric elements m and f3 , can be determined from eq. 
([?]) and (JsJ) , up to some functions of the time-like coordinate u. In the context 
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of this approach the radial dependence of p and P is assumed to be the same 
as in the static "seed" solution. 

The rationale behind the assumption on the r dependence of the effective 
variables P and p, can be grasped in terms of the characteristic times for different 
processes involved in a collapse scenario. If the hydrostatic time scale Thydr, 
which is of the order ~ \J\JUp (where G is the gravitational constant and 
p denotes the mean density) is much smaller than the Kelvin- Helmholtz time 
scale ( Tkh ), then in a first approximation the inertial terms in the equation 
of motion can be ignored Jl^| . Therefore in this first approximation (quasi- 
stationary approximation) the r dependence of P and p are the same as in the 
static solution. Then the assumption that the effective variables ^) and (^|) 
have the same r dependence as the physical variables of the static situation, 
represents a correction to that approximation, and is expected to yield good 
results whenever Tkh 3* Thydr- Fortunately enough, Tkh Thydr, for 
almost all kind of stellar objects. 

Those functions of the time-like coordinate u that remain arbitrary can be 
obtained from a system of ordinary differential equations ( The System of Surface 
Equations) emerging from the field equations evaluated at the boundary surface 
and the junction conditions. 

The equations corresponding to the junction conditions (eq. (49), and (50) 
in @) lead to 

Pxa (1 + o) = /(u)a (9) 

and 

fhia (1 + a) = g(u)a . (10) 

Where f(u) and g(u) are arbitrary functions of their arguments, differentiation 
with respect to u and r are denoted by subscripts and 1, respectively and the 
subscript a indicates that the quantity is evaluated at the boundary surface. 
Finally, in radiation coordinates, d takes the form 



a — 



da u xa A 2m a 



du 1 — U); 

Also, from the junctions conditions (eq. (40) in [|l2|) we get 



1-— . (11) 



or equivalently 



2a/3 lQ ( 1 + d-— ) =m lQ (13) 



From equations (|13[) , d9j) and (|10|) we obtain an expression relating f(u) and 
g(u), namely 

2a(l + a-^jf(u)=g(u). (14) 



G 



However, it can be checked by simple inspection that neither the field equations 
nor the junction conditions impose further restrictions on these functions of u, 
therefore one of them remains completely arbitrary for each model. 
Expanding (||) for a <§; 1 , we get 

1 1 1 

e = —a- - — 

r, 2 r J 



n - -~a 3 H , (15) 



thus, as expected, up to first order the eccentricity is proportional to a. Now 
equations (|l(]) and (|l4|) lead to 



miq (1 + a) 
2a (l + a-^)/( U ) 



a - ^ TZTA (16) 



or, using eq. (0), 



a = 2 .; ° 2^T";.. • ( 17 ) 



nap a (u) (1 + d) 
(l + 4 _2^)/(tt) 

It is easy to see that two different types of variables arc needed to determine 
a from (|l7j). On one hand, we notice a,d,/5 a , and rh a which can be obtained 
from the integration of the Surface Equations and the particular "seed" solution 
selected to be modeled. On the other hand, as we have stated before, f(u) 
remains completely arbitrary and its choice completes the characterization of 
the model. 

In the next section we shall present three different families of models obtained 
from different static equations of state. 



3 The Modeling 

In order to show that there is no model independent tendency for the eccentricity 
to attain a maximum, we shall consider three different families of equations of 
state. Each one will be modeled with three different choices of the arbitrary 
function f(u). We shall work out three models, previously studied for the 
spherical (nonrotating) case: Schwarzs child- like Jl4|, 15 ,Tolman-VI-like p"5| 



and Tolman-V-like |14| , |16| , |17| solutions. In the static limit the Schwarzschild- 
like homogeneous solution represents an incompressible fluid of constant density. 
Static Tolman VI solution approaches that of a highly relativistic Fermi Gas 
and, therefore, with the corresponding adiabatic exponent of 4/3. And, for the 
static Tolman V solution, the radiation relation P/p ~ 1/3, is maintained at 
the center of the distribution during the contraction. 

All the above equations of state will be used as a "seed" solution to study 
the evolution of the eccentricity. The first family to be considered is the 
Schwarzschild solution. It is the same example presented in ref. |p^| . 

3 (1-F(u)) 

P = k{u) = ^ A{u)* ' (18) 
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and 



3(1- F(u))(l + F(u) (fi(u)-l)) 



(19) 



AA(u) F(u) Q(u) f(u) 

where the exterior radius a, the total mass m, and the timelike coordinate 
u are scaled by the initial total mass m(u = 0) = m(Q), i.e. A = a/m(0); 
M = m/m(0) and u — u/m(Q). We have also defined 



and 



F(u) 



i 



2M(u) 



1 



(20) 



(21) 



1 — LOxa 

The second case is the slowly rotating Tolman-V-like models. In this case we 



have 



where 



with 



and 



P 



1 (8{u) 



z{u)r~~ 



8vr V r 2 

{l-F{u)){l + F{u) (Sl(u)-1)) 
3 A(u) F(u) /(«) 



6(u) = ~{F(u)-l)(2S}(u)-5) 



z{u) = 
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(l-F(u))(4n(u)-3). 



(22) 
(23) 
(24) 
(25) 



21 A{u) 7 ^ 

Finally we take the Tolman Vl-like family of models, with the effective den- 
sity represented by 



P 



3h(u) 1 (1- F(u)) 
r 2 8ir r 2 



and for the eccentricity we have found 

(l-F(u))(l + F(u) (Q(u)-1)) 
4A(u) F(u) Q(u) f(u) 



(26) 



(27) 



As we have mentioned above, the evolution of these three families of eccen- 
tricities (|i"9|), (p7|), ( p3| ) have been calculated for three different expressions of 
f(u), namely: 

cte = 1 
1 +u 



/(«) 



l + Cexp-(^^)" 



(28) 



with Utrig a triggering parameter, C a modulation constant and a the width of 
the pulse. 
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Simulations have been carried out for a large number of sets of initial data 
and emission pulses. We shall present in the next section only the most relevant 
cases. 

The radiation flux has been assumed to be in the free streaming out limit 
approximation, i.e. neutrinos and/or photons mean free path is of the order of 
the dimensions of the spheroid. It is worth mentioning that for the modeling 
performed, the shape and the intensity of the emission pulse does not seem 
to play an important role, we shall only exhibit the results obtained with a 
Gaussian pulse centered at u = u p 

■ AM rad 1 (u-U p \ 2 

with A the width of the pulse and AM ra ^ the total mass lost in the process. All 
models presented in the next section have been simulated using 

AM rad = 0.001 A = 0.1 u p = 15 (30) 

The sets of initial data relevant to be discussed are: 



Schwarzschild-like (displayed in Figures 1 and 2) 
M(0) 



A(0) 

and 

M(0) 



0.185; fi(0) = 0.933; A(0) = 5.404; F(0) = 0.630 (31) 



= 0.303; fi(0) = 0.933; A(0) = 3.300; F(0) = 0.394 (32) 



A(0) 

• Tolman V-like (displayed in Figure 3 ) 
M(0) 



0.217; fi(0) = 0.809; A(0) = 4.600; F(0) = 0.565 (33) 



A(0) 

• Tolman Vl-like (displayed in Figure 4 ) 
M(0) 



A(0) 



0.270; n(0) = 1.000; A(0) = 3.700; F(0) = 0.459 (34) 



The running time for the models (values for u) are controlled by the physical 
relevance of the variables involved. As we shall see in the next section the 
behavior of a is highly model dependent, no matter how compact the object 
becomes. 
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4 Discussion of Results 



Let us start by analyzing the Schwarzschild-like model. Figures la through 
lc display the evolution of the eccentricity for the above mentioned choices of 
f(u) in equation (jH) with the set of initial data (|3l|). Figure Id indicates the 
evolution of F(u). Observe that cases 



= { l + Cexp-(^-) 2 < 35 > 

yield results very similar to those reported in although we are considering only 
first order terms ! These result are presented in Figures la and lc, respectively. 
From Figures la, lc and Id we can see that the eccentricity for the above two 
cases present a maximum when the surface gravitational potential reaches 

3$ -|p -«-»-".. m 

The other model, corresponding to the case /(it) = 1 + u also gives a stationary 
value for the eccentricity for the above value for the gravitational potential at 
the surface of the configuration. 

Thus, up to the first order this model would seem to confirm the result 
found in @, ||. However, this is not so. Indeed, Figures 2a to 2d represent the 
evolution of the eccentricity, e(it) , and F(u) , for another set of initial data (see 
equations (|32|)). Now the stationary value for e(u) is obtained for ^jf-f ~ 0.33 



. It is worth noticing that for a given set of initial data the final value of e(it) 
is independent of the f(u) . 

Let us now turn to consider the Tolman V-like family of models with the 
set of initial data (|33|). Figures 3a through 3c exhibit the evolution of the 



eccentricity concerning the functions /(it) in equation (28). In Figure 3d the 
evolution of F(u) is displayed. As it is apparent from these figures, there is no 
maximum in the value of the eccentricity, even though the surface gravitational 
potential goes well beyond the supposed critical value of ^jj^y ~ 0.18 (see Figure 
3d). 

Finally, Figures 4a to 4d refer to Tolman VFlike family of models, and point 
in the same direction as the previous case. In fact, Figures 4a to 4c show that 
either the eccentricity grows up continuously (/(it) is one of the above eq. (pq)) 
or attains a minimum and grows up again (f{u) = 1 + u), even though the 
surface gravitational potential may reach values as high as ~ 0.31 (Figure 

4d). 

Summarizing our results. It is clear that, up to the first order in angular 
parameters, the evolution of the eccentricity is highly model dependent. If 
second order contributions are going to be considered, they should be assumed 
much smaller than the first order ones. Therefore, it is not reasonable to expect 
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that second order contributions could bend the curves displayed in Figures 3a, 
3c, 4a and 4c as to create a maximum for e(u) or that they can subdue the 
dependence of the final value of e(u) on the initial data for the Schwarzschild- 
like models. In other words, we can not see how the second order effects might 
affect qualitatively the model dependence exhibited here. 

Of course, it remains to be seen whether the effects described in & ||] appear 
in a model independent way in the exact theory of relativistic rotating bodies. 
For that purpose, however, we would need an exact interior solutions for the 
axially symmetric source, which, as it is well known, is not available yet. 

We would like to conclude with the following comment. Although the main 
purpose here has been to bring out the model dependence of the evolution of the 
eccentricity, and not specific modeling of rotating sources, it should be noticed 
that in the examples above the energy density is always positive and larger 
than pressure everywhere within the fluid distribution and the fluid velocity as 
measured by the locally Minkowskian observer is always less than one. Although 
we do not include figures of physical variables (some of them for the case / = 1 
may be found in reference |L2|), the reader may convince himself of that by 
noticing that the three models considered here have been studied in detail in 
the spherically symmetric case [H| [l7| , in which all models have a correct 
physical behaviour. Of course we do not expect that first order perturbation 
may change that situation. 

Two of us (H.H. and L.A.N) gratefully acknowledge the financial support by 
the Consejo de Desarrollo Cientffico Humam'stico y Tecnologico de la Universi- 
dad de Los Andes, under project C-720-95-B and by the Programa de Formacion 
e Intercambio Cienti'fico (Plan II). 
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Figure Captions 



Figure 1 The evolution of the eccentricity for the Schwarzschild-like model 
with ^fa) — 0.185. Figures la through lc correspond to functions: 

f(u) = cte = l;f(u) = 1 + m: and /(u) = 1 + Cexp - (^^^) 2 , re- 
spectively. Figure Id indicates the evolution of the gravitational potential 
at the surface: F(u) = 1 — 2 ^ff • 

Figure 2 The evolution of the eccentricity for the Schwarzschild-like model 
with = 0.303. Figures 2a through 2c correspond to functions: f(u) = 

f(u) = 1 + it; and /(it) = 1 + C exp — ^ "~^ tr, g ^ , respectively. Figure 
2d indicates the evolution of the gravitational potential at the surface: 

Figure 3 The evolution of the eccentricity for the Tolman V-like model with 
= 0.217. Figures 3a through 3c correspond to functions: /(it) = cte = 

f(u) = 1 + it; and /(it) = 1 + C exp — ^ - u * r13 ^ , respectively. Figure 
3d indicates the evolution of the gravitational potential at the surface: 
*•(«) = !- ^ 

Figure 4 The evolution of the eccentricity for the Tolman Vl-like model with 
= 0.270. Figures 4a through 4c correspond to functions: /(it) = cte = 

/(li) = 1 + it; and /(it) = 1 + C exp — ^ "~^ tr, g ^ , respectively. Figure 
4d indicates the evolution of the gravitational potential at the surface: 



13 



Eccentricity 
0.0151 



40 60 
u 



Eccentricity 
0.014 




40 60 
u 



(b) 



Eccentricity 
0.0151 




20 



40 60 
u 



40 60 
u 



(d) 



Figure 1: 



